MAPYCUS MAXIMUS

grid <- create_test_grid(range = c(-1, 1), spacing = 0.1)
transform <- fisheye_fgc(grid, r_in = 0.3, r_out = 0.5, zoom_factor = 1.5, squeeze_factor = 0.5, method = "outward")
plot_fisheye_fgc(grid, transform, r_in = 0.3, r_out = 0.5)

grid_df <- as_tibble(grid)
transform_df <- as_tibble(transform) |>
  dplyr::mutate(
    zone = attr(transform, "zones"),
    r_orig = attr(transform, "original_radius"),
    r_new  = attr(transform, "new_radius")
  )
arrows_df <- cbind(grid_df, transform_df)

ggplot() +
  # Draw arrows showing movement
  # Original points
  geom_point(
    data = grid_df, aes(x = x, y = y),
    size = 0.6, alpha = 0.7, color = "black"
  ) +
  # Transformed points
  geom_point(
    data = transform_df, aes(x = x_new, y = y_new, color = zone),
    size = 0.6, alpha = 0.7
  ) +
    geom_segment(
    data = arrows_df |> filter(zone != "context"),
    aes(x = x, y = y, xend = x_new, yend = y_new, color = zone),
    arrow = arrow(length = unit(0.02, "npc")),  # optional arrowhead
    alpha = 0.6, size = 0.5
  ) +
  scale_color_manual(values = c("focus" = "#c60000ff", 
     "glue" = "#141497ff", 
     "context" = "#FFCC00")) +
  coord_equal() +
  theme_minimal(base_size = 14) +
  labs(title = "Fisheye Transformation: Point Movement", x = "x", y = "y")

library(purrr)
library(dplyr)
library(stringr)

zoom_seq <- seq(1, 20, by = 0.1)
center_pt_proj <- melbourne

# Calculate melbourne centroid once
melbourne_cent <- melbourne |> st_union() |> st_centroid()

fisheye_frames <- map_dfr(zoom_seq, function(z) {
    # Apply fisheye transformations
    vic_fish_new  <- sf_fisheye(vic_fish, center = center_pt_proj,
                          r_in = 0.34, r_out = 0.5, zoom_factor = z)
    conn_fish <- sf_fisheye(conn_small, center = center_pt_proj,
                          r_in = 1.12, r_out = 5, zoom_factor = z)
    hosp_fish <- sf_fisheye(hosp_points, center = center_pt_proj,
                          r_in = 1.23, r_out = 2.3, zoom_factor = z)
    racf_fish <- sf_fisheye(racf_points, center = center_pt_proj,
                          r_in = 1.58, r_out = 2.3, zoom_factor = z)
    
    # Transform melbourne_cent to match coordinate system
    melbourne_cent_trans <- melbourne_cent |> st_transform(st_crs(hosp_fish))
    
    # Calculate distances
    distances_racf <- st_distance(racf_fish, melbourne_cent_trans) |> as.numeric() |> round()
    distances_hosp <- st_distance(hosp_fish, melbourne_cent_trans) |> as.numeric() |> round()
    
    # Add distances and counts
    racf_fish <- racf_fish |> mutate(dist = distances_racf) |> add_count(dist)
    hosp_fish <- hosp_fish |> mutate(dist = distances_hosp) |> add_count(dist)
    
    # Filter based on distance frequency
    racf_fish_filter <- racf_fish |> filter(n != max(n))
    hosp_fish_filter <- hosp_fish |> filter(n != max(n))
    
    # Filter connections
    conn_fish_filter <- conn_fish |> filter(source %in% racf_fish_filter$source)
    conn_fish_filter <- conn_fish_filter |> arrange(desc(weight)) |> head(10)
    
    # Filter based on connections
    racf_fish_filter <- racf_fish_filter |> filter(source %in% conn_fish_filter$source)
    hosp_fish_filter <- hosp_fish_filter |> filter(destination %in% conn_fish_filter$destination)
    
    # Shorten hospital names
    hosp_fish_filter <- hosp_fish_filter %>%
      mutate(hosp_name_short = hosp_name %>%
        str_remove(" Health Service$") %>%
        str_remove(" Hospital$") %>%
        str_remove(" Inc\\.$") %>%
        str_remove(" Ltd\\.$") %>%
        str_remove(" Private$") %>%
        str_remove(" Centre$") %>%
        str_remove(" Campus$") %>%
        str_remove("^The ") %>%
        str_replace(" Private Hospital", "") %>%
        str_replace(" Health$", "") %>%
        str_replace(" Rehabilitation Hospital", " Rehab") %>%
        str_replace(" Rehabilitation Centre", " Rehab") %>%
        str_replace(" Day Surgery", " Day Surg") %>%
        str_replace(" District Hospital", "") %>%
        str_replace(" Regional Hospital", " Regional") %>%
        str_squish()
      )
  
    tibble(
      zoom_factor = z,
      vic = list(vic_fish_new),
      conn = list(conn_fish_filter),
      hosp = list(hosp_fish_filter),
      racf = list(racf_fish_filter)
    ) 
})

# Reshape to long format
fish_long <- map_dfr(1:nrow(fisheye_frames), function(i) {
  z <- fisheye_frames$zoom_factor[i]
  
  bind_rows(
    fisheye_frames$vic[[i]]  %>% mutate(type = "vic",  zoom_factor = z),
    fisheye_frames$conn[[i]] %>% mutate(type = "conn", zoom_factor = z),
    fisheye_frames$hosp[[i]] %>% mutate(type = "hosp", zoom_factor = z),
    fisheye_frames$racf[[i]] %>% mutate(type = "racf", zoom_factor = z)
  )
})

library(gganimate)
library(ggplot2)

p <- ggplot() +
    # VIC map Polygon
    geom_sf(data = subset(fish_long, type == "vic"),
        fill = NA, color = "grey") +
    # Hospital labels (using short names)
    geom_sf_label(data = subset(fish_long, type == "hosp"),
        aes(label = hosp_name_short), color = "red", size = 2.5, alpha = 0.5) +
    # RACF labels
    geom_sf_label(data = subset(fish_long, type == "racf"),
        aes(label = racf_name), color = "blue", size = 2.5, alpha = 0.5) +
    # Hospitals (points)
    geom_sf(data = subset(fish_long, type == "hosp"),
        color = "red", size = 2) +
    # Age Care Facilities (points)
    geom_sf(data = subset(fish_long, type == "racf"),
        color = "blue", size = 2) +
    # Connections
    geom_sf(data = subset(fish_long, type == "conn"),
        aes(alpha = weight), color = "black") +
    coord_sf(crs = st_crs(fish_long)) +
    labs(title = "Transportation between Hospital and Age Care Facilities in Victoria during COVID-19",
         subtitle = "Zoom: {current_frame}×") +
    theme_map() +
    theme(
      plot.title = element_text(size = 15, hjust = 0.5, margin = margin(t = 10, b = 5)),
      plot.subtitle = element_text(size = 10, hjust = 0.5),
      plot.margin = margin(t = 20, r = 10, b = 10, l = 10)
    ) +
    transition_manual(zoom_factor)

anim <- animate(
  p,
  fps = 25,
  duration = 8,
  width = 1366,
  height = 768,
  res = 150
)

anim_save("fisheye_zoom_gganimate.gif", animation = anim)